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, Abstract 
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Homoclinic snaking refers to the sinusoidal "snaking" continuation curve of homoclinic orbits near a 
heteroclinic cycle connecting an equilibrium E and a periodic orbit P. Along this curve the homoclinic 
orbit performs more and more windings about the periodic orbit. Typically this behaviour appears in 
reversible Hamiltonian systems. Here we discuss this phenomenon in systems without any particular 
structure. We give a rigorous analytical verification of homoclinic snaking under certain assumptions on 
the behaviour of the stable and unstable manifolds of E and P. We show how the snaking behaviour 
depends on the signs of the Floquet multipliers of P. Further we present a nonsnaking scenario. Finally 
we show numerically that these assumptions are fulfilled in a model equation. 
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1 Introduction 

In this paper we study analytically a certain continuation scenario, the so-called Homoclinic Snaking, of 
homoclinic orbits in systems without particular structure such as reversibility or Hamiltonian structure. 
, Consider an ordinary differential equation x = f(x) with x G W l . Given an equilibrium E, a homoclinic 
orbit to E is a solution that converges to E as t — > ±oo. If E is hyperbolic, a homoclinic orbit lies in 
the intersection of the stable and unstable manifolds W S (E) and W U (E) of the equilibrium. For general 
differential equations, stable and unstable manifolds of hyperbolic equilibria will typically not intersect by 
the Kupka-Smale theorem [22]. Transversality arguments show that in one-parameter families of differential 
equations one can expect an intersection, and hence a homoclinic orbit, to occur persistently at an isolated 
parameter value. Therefore a continuation of a homoclinic orbit can be carried out in two-parameter families 
of differential equations. 

Indeed, the homoclinic orbits we consider live within a small neighbourhood of a heteroclinic cycle connecting 
a hyperbolic equilibrium E and a hyperbolic periodic orbit P. Such a cycle consists, besides E and P. of 
two orbits 7 EtoP and 7 PtoE with lini{_>._ o0 7 Eto p(i) = E, lim^oo 7 EtoP (t) = P, and lim^-oo 7 PtoE (t) = P, 
limt-^oo 7p toE (t) = E, respectively The orbit 7 EtoP is called heteroclinic orbit connecting E to P or in short, 
EtoP connecting orbit , or just EtoP connection. A similar terminology holds for 7 PtoE just with interchanging 
E and P. The complete cycle we call EtoP cycle. More precisely the homoclinic orbits to the equilibrium 
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under consideration are one-homoclinic orbits with respect to the given EtoP cycle. That means they move 
once along the cycle before returning to the equilibrium. 

Replace the periodic orbit P by a hyperbolic equilibrium E, and assume that for a critical value in a 
two-dimensional parameter space there exists a heteroclinic EtoE cycle. Then, typically there is a curve 
emanating from the critical value such that for all parameter values on this curve there exits a homoclinic 
orbit to E [7, 27]. The homoclinic orbit spends more and more time near E when moving along this curve 
towards the critical point. 

Considering homoclinic orbits to E in a neighbourhood of an EtoP cycle, we also find that along the 
continuation curve the homoclinic orbits spend more and more time near the periodic orbit. In contrast to 
the EtoE-case the continuation curve does not converge to a point, but it approaches a curve segment as h\ 
does in Figure 2 below. This makes the consideration global in the parameter space. Before discussing this 
in more detail we consider the problem from the homoclinic snaking point of view. 

In the context of ordinary differential equations the notion Homoclinic Snaking originally denotes a continu- 
ation scenario of homoclinic orbits in reversible Hamiltonian systems. In Hamiltonian systems the situation 
is somewhat different to the one described above. Both the stable and the unstable manifold of a hyperbolic 
equilibrium are in the same levelset of the Hamiltonian. Therefore they will typically intersect transversely 
(within this levelset) and a homoclinic orbit can be expected to occur persistently in a single system. Hence 
a continuation can be done in one-parameter families. Typical continuation curves related to a homoclinic 
snaking scenario are displayed in Figure 1, cf. [1, Figure 1.1]. 




0.64 0.7 ft 

Figure 1: Bifurcation diagram for homoclinic orbits of the steady states —(1 + d^) 2 U — fiU + 2U 3 — U 5 = of the 
Swift-Hohenberg equation. The central panel shows the typical snakes and ladder structure. The outer panels show 
the U vs. time plot for the corresponding points on the snaking curves. 

Figure 1 shows continuation curves of homoclinic orbits related to the fourth order scalar equation 

- (1 + dl) 2 U- U.U + 2U 3 - U 5 = 0. (1.1) 

Considered as a system in K 4 , equation (1.1) is a reversible Hamiltonian system with family parameter //. 
The two intertwined wiggly curves correspond to symmetric (w.r.t. the reflection x i— > —x and U i— > —U) 
homoclinic orbits asymptotic to a hyperbolic equilibrium. These curves are also called snaking curves. The 
term snaking is due to the sinusoidal shape of the continuation curves. Indeed, the homoclinic orbits under 
consideration are one-homoclinic w.r.t. a symmetric heteroclinic cycle connecting a symmetric equilibrium 
with a symmetric periodic orbit. Restricted to the corresponding levelset of the Hamiltonian, the periodic 
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orbit is hyperbolic. The excursion of the homoclinic orbit to the periodic orbit lasts longer and longer along 
the continuation curves. The homoclinic orbit performs more and more windings about the periodic orbit 
- this corresponds to the increase of its L 2 -norm. Roughly speaking, the /x-range of the snaking curves is 
the /x-range for which the heteroclinic cycle does exist. At the cndpoints of the /x-interval the involved EtoP 
and PtoE connections (note that they are images of each other by the reversing symmetry) simultaneously 
undergo a saddle-node bifurcation. Note that the snaking curves indicate saddle-node bifurcations of the 
symmetric homoclinic orbits. These bifurcation points are close to the endpoints of the ix-range for the 
heteroclinic cycle. In addition to the snaking curves, Figure 1 also displays a ladder structure. The rungs 
connecting the two snaking curves correspond to asymmetric homoclinic orbits to the equilibrium. The 
asymmetric homoclinic orbits bifurcate from the symmetric ones via pitchfork bifurcation. These bifurcation 
points are close to the saddle-nodes of the symmetric homoclinic orbits. However, in the context of the present 
treatment it is enough to focus on the features of one single snaking curve. Bifurcation diagrams as displayed 
in Figure 1 have been discussed for instance in [1, 3, 9, 28]. For a more complete list of references we refer 
to [1], but with the addition that homoclinic snaking, also called collapsed snaking, near EtoE cycles in 
one-parameter families of reversible systems has been studied in [16, 17]. 

More recently, Krauskopf, Oldeman and Ricfl [19, 20] numerically discovered a similar effect in a system 
without any particular structure such as reversibility or Hamiltonian structure. The corresponding family 
of vector fields in R 3 with family parameter {v\, v^) can be written in the form 



x = v\x — y + x sin Lp — (x 2 + y 2 )x + 0.01(2 cos ip + v^) 2 
i] = v\y + x + y sin if — (x 2 + y 2 )y + 0.0l7r(2 cos ip + v^) 2 
(p= v-i — (x 2 + y 2 ) + 2cos<p 



F(x,y, ip,v u v 2 ). 



(1.2) 



Figure 2 displays a continuation curve h\ for a homoclinic orbit detected numerically in [19, 20]. Again, and 
not only due to its shape, we address this curve as snaking curve. As in the reversible case, the homoclinic 
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Figure 2: Snaking diagram of a three-dimensional laser model (1.2). Panel (a) shows the snaking curve h\ together 
with the continuation curve Cb of the PtoE connecting orbits and the locus tb of saddle-nodes of the EtoP connecting 
orbits. Panels (b) - (d) show y vs. time plots of the y component at the corresponding points (b) - (d) in panel (a). 



orbit under consideration is one-homoclinic w.r.t. a heteroclinic cycle connecting a hyperbolic equilibrium 
with a hyperbolic periodic orbit, and along the continuation curve the homoclinic orbits performs more 
and more windings about the periodic orbit - cf. panels (b)-(d) in Figure 2, which show plots of one state 
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variable corresponding to points indicated in panel (a). The plot of the L 2 -norm of the (x,?/)-part of the 
solution versus the parameter v\ or ^2, respectively, behaves as in the Hamiltonian case, cf. Figure 3. The 
snaking behaviour w.r.t. both parameters is due to the declination of the curve q, in Figure 2(a). This 
defines intervals within which the parameters move while h\ approaches q,. We refer also to Section 6 for 
more numerical results regarding this system. 
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Figure 3: Snaking diagrams of a three-dimensional model (1.2). Shown are plots of the I? -norm of (x,y) vs. v\ and 
V2, respectively, along the snaking curve h\. 



Although system (1.2) has the same dimension as the restriction of the above Hamiltonian system to a 
levelset, the geometry is quite different. Here one of the heteroclinic connections constituting the EtoP cycle 
does not lie in a transversal intersection of the corresponding stable and unstable manifolds. Hence, by the 
same arguments as given above for homoclinic orbits, one can expect to find it on a continuation curve in 
the i^-space. This is the curve q, in Figure 2(a). The other connecting orbit is robust - as the ones 
in the Hamiltonian case. This connecting orbit exists within the stripe delimited by the curves tb, cf. again 
Figure 2(a). In other words, the region where the snaking curve is located is related to the existence of 
the robust heteroclinic connection between the equilibrium and the periodic orbit. Remarkably, the snaking 
curve accumulates at the curve segment defined by the intersection of the curve Cb with the stripe delimited 
by the curves tb- In other words, the snaking curve accumulates at the line segment for which the EtoP cycle 
does exist. 

In this paper we give an analytical verification of the bifurcation diagram in Figure 2 within a more general 
setup. We consider a general two- parameter family of differential equations in R . In this general context we 
denote the family parameter by (Ai, A2) taken from a closed rectangle Q. We assume that there is a closed 
interval I2 such that for all (Ai,A2) G {0} x I 2 there exists an EtoP cycle connecting an equilibrium and 
a periodic orbit, both are assumed to be hyperbolic, cf. Hypothesis 2.1. Let W S (E) be two-dimensional. 
We want to note that this setting is related to the vector field — F in (1.2). Compare also the caption of 
Figure 2. Further we constitute conditions concerning the global intersection behaviour of W U (P), W S (E) 
and W U (E), W S (P), respectively, cf. Hypotheses 2.3 - 2.6. Note that all manifolds depend on A, which is 
so far suppressed from the notation. Indeed, the snaking behaviour is mainly influenced by the behaviour of 
the intersection of W U (P) and W S (E). Consider a small torus T around P, and let E ou * be a small stripe 
on this torus around W U (P) H T ■ Similarly we define £*" as a small stripe on this torus around W S (P). In 
r C S 1 x Q we collect the intersections of W U {P) and W S (E) in E ou * depending on A: 

r := {(if, A) : W S (E, A) n iy M1 (P(</?, A), A) n S out ^ 0}, 

where W uu (P((p, A)) is the strong unstable fibre of P(tp,\) 6 P. The assumption that T is graph of a 
function z = z(ip, Ai) is essential for the snaking behaviour. Our hypotheses on z imply amongst others that 
at the endpoints of I2 the PtoE connection undergoes a saddle-node bifurcation, and what is more, for each 
A G {0} x int I2 there are at least two EtoP cycles. 

To get a better idea of what these assumptions include, suppose for simplicity that the trace of W S (E) in 
S otl * is a closed curve which is simple over W U {P). So W S (E) can be seen as function Z of values in W U {P). 
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To simplify matters further we assume that by changing A these curves will be shifted against each other 
without changing their shape. Then Z and z(-,0) are directly related. 

For Ai = and all A2 G I21 there is a unique EtoP connection; all those A2 can be written in the form 
A2 = z((p, 0). Our hypotheses on the EtoP connection provide the existence of a 27r-periodic function (Pq{L) 
defining the base point of the strong stable fibre of P which intersects in E m the EtoP connection related 
to Ai =0, A 2 =z(<p%(L) + 2L,0). 

Our main snaking result, cf. Theorem 2.1, says that all one-homoclinic orbits near the primary EtoP cycles 
lie on one continuation curve, which we refer to as the snaking curve. This curve can be parametrised by 
the flight time L of the one-homoclinic orbit between E m and T, out . For L — > 00, this curve accumulates at 
{0} x li- It turns out that the shape of the snaking curve is mainly determined by z - more precisely \i(V) 
arises as a perturbation of graph z((p^(L) + 2L, 0). 

In the scenario covered by Theorem 2.1 the nontrivial Floquet multipliers of P are positive. Indeed, the 
snaking behaviour depends on the sign of these multipliers. If they are positive, the local (un)stablc manifold 
of P is topologically a cylinder, while for negative multipliers these local manifolds are topologically a Mobius 
strip. Therefore, for positive multipliers both E m and consist of two connected components. In our 
analysis however, only one of these components, in each case, plays a role - and there is only one way for the 
transition from S m to T, out . But if the multipliers are negative, then E m and S°" 4 are connected, and both 
are winding twice around P on T. This results in two different ways for the transition from E m to S out , and 
this causes the existence of two snaking curves approaching {0} x I2 from different sides, cf. Theorem 4.1. 

Next we abandon our assumption on T being graph of a function and assume rather that for fixed Ai the 
set r is a closed curve, and replace Hypothesis 2.4 by Hypothesis 5.1. This prevents snaking. In this case 
no longer all homoclinic orbits are on one continuation curve, instead there exists a sequence of closed 
homoclinic continuation curves in the A-plane accumulating at {0} x I2, cf. Theorem 5.1. 

Numerically we show that for the motivating system (1.2) our hypotheses generating snaking are fulfilled. 
In particular we verify Hypothesis 2.4. 

For our analysis we use Fenichel coordinates near the periodic orbit, and within this setting we consider 
solutions of a Shilnikov problem which we glue together with the stable and unstable manifolds of the 
equilibrium. This procedure is the same as the one utilised in [1] and [14] to study reversible Hamiltonian 
systems. To our knowledge there are only a few further works presenting analytical results for the dynamics 
near EtoP cycles. In [15] a Lin's method approach has been extended to treat heteroclinic chains involving 
period orbits. These results are applied to EtoP cycles, in particular to detect nearby one-homoclinic orbits. 
However, the results are more local in nature. These results concern, in the context of Figure 2, the existence 
of one-homoclinic orbits for parameter values in the neighbourhood of certain points on the curve Cf, - but 
not in the neighbourhood of an entire segment of Cb as in the present paper. In [23, 24] a somewhat different 
(in handling the flow near P) Lin's method approach has been used to study EtoP cycles of codimension-one 
and codimension-two. In this language the EtoP cycle considered in the present paper are of codimension- 
one. Rademacher's results about homoclinic orbits near codimension-one EtoP cycles are of the same nature 
as the ones in [15]. 

In our analysis we assume the existence of a primary EtoP cycle, and we make assumptions about its global 
(in parameter space) behaviour. In the Swift-Hohenberg equation the existence of heteroclinic connections 
has been investigated analytically in [6, 18]. 

In [4], amongst others, homoclinic snaking caused by an EtoP cycle in systems in R 3 is considered. Using 
a combination of geometric and analytical arguments, the snaking behaviour as displayed in Figure 2 is 
explained. More precisely, based on a leading term approximation of the bifurcation equation, one-homoclinic 
orbits near the saddle-node points of the EtoP cycles are determined. Then the entire snaking curve is 
deduced by using geometric arguments. This has been done for both, when P has positive or negative 
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Floquet multipliers. Here in the present paper on the contrary, we give a rigorous analytical verification of 
these scenarios. 

Numerically the homoclinic snaking scenario in the addressed model, which is inspired by semiconductor laser 
dynamics, has been considered in several papers. The snaking curve h\, cf. Figure 2, was first revealed in [19]. 
In [20] the relation to the organising EtoP cycle was numerically discovered in a bifurcation diagram similar 
to that in Figure 2. This system was further investigated in [4, 15]. In [4, 5] a similar snaking behaviour 
was (numerically) observed in a nine-dimensional model equation of intracellular calcium dynamics. The 
remarkable feature in the bifurcation diagram is that turning points of the snaking curve accumulate on six 
different values (and not on only two as in the one displayed in Figure 2). Note that our analysis is carried 
out only for the case of a three-dimensional state space. However, the set T remains a curve also in higher 
dimensions. Then a corresponding function z can be defined, and the addressed feature can be explained by 
the number of critical points of z, cf. also Figure 9 below. 

In [2] homoclinic snaking in the transition from reversible Hamiltonian systems to general systems using 
the example of the Swift-Hohenbcrg equation is considered. Indeed this happens in R , but the observed 
snaking or nonsnaking behaviour in the perturbed system, respectively, discloses features we discuss for 
general systems in R 3 . In [2, Figure 6] snaking curves of two different homoclinic orbits are shown - each 
curve displaying a behaviour as shown in Figure 3 of the present paper. These two different homoclinic 
orbits are remains of the unperturbed reversible Hamiltonian system. Isolas of homoclinic orbits as shown 
in [2, Figures 2 and 3] are discussed in Section 5 of this paper. In symmetric systems those isolas are also 
observed in [3, Figure 24]. However, it is not the aim of this paper to explain those transition processes. 

This paper is organised as follows. In Section 2 we present our hypotheses and formulate the main snaking 
result, Theorem 2.1. The proof of Theorem 2.1, is then carried out in Section 3. In Section 4 we treat 
negative Floquet multipliers. The results are summarised in Theorem 4.1. Afterwards we discuss one 
possible nonsnaking scenario in Section 5; Theorem 5.1 covers the results of this section. In Section 6 we 
verify numerically Hypothesis 2.4, the main snaking assumption, in the laser model (1.2). 

2 Setup and main results 

We consider a smooth family of differential equations 

x = f(x,\), xeR 3 , A = (Ai, A 2 ) e Q C K 2 , (2.1) 

where Q = J\ x J 2 is a closed rectangle with 6 intQ; J\, J-2 are closed intervals. 
We assume the following 

Hypothesis 2.1. (i) /(0,A) = 0; The equilibrium E := {x = 0} is hyperbolic, and dimW u (E, A) = 1, 
dimW s (.E,A) = 2. 

(ii) For all A 6 Q there is a hyperbolic periodic orbit P. Further let dimW u (P, A) = 2, dimW s (P, A) = 2. 
For all A the minimal period of P is 2ir. 

(Hi) There is a maximal interval I2 '■= [A 2 1-^2] $! J2, A2 < ^2, such that for A <E {0} x I2 there is a 
heteroclinic cycle connecting E and P. 

The constant minimal period can always be achieved by an appropriate time transformation. The interval 
I2 is maximal in the sense that for (Ai = 0, A2) and A2 > A2 or A2 < A2 there in no complete cycle. 
More precisely with our choice of dimensions, typically the EtoP connection is of codimension-one - that 
means it appears along a curve in parameter space. This curve is the A2-axis and the connection splits up 
when moving off the A2-axis. On the other hand, the PtoE connection is typically robust. Nevertheless, 
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by changing parameters within a wider range this connection can disappear, for instance in the course of a 
saddle-node bifurcation. These scenarios are made more precise by additional hypotheses below. 

The three-dimensional state space enforces that both nontrivial Floquet multipliers of P have the same sign, 
cf. [26]. 

Hypothesis 2.2. The nontrivial Floquet multipliers of P are positive. 

The positivity of the Floquet multipliers of P is exploited in the Fenichel normal form near P, cf. Lemma 2.2 
below. However, in Section 4 we relaxe this hypothesis. 

The following lemma can be seen as a motivation for our further considerations. Roughly speaking, it says 
that under certain transversality conditions on each curve k intersecting {Ai = 0} transversely, there is a 
sequence of parameter values accumulating at {Ai = 0} for which a onc-homoclinic orbit to the equilibrium 
does exist, cf. [15, Corollary 4.3]. 

Lemma 2.1 ([15]). Assume Hypotheses 2.1 and 2.2, and let k = k(/i) be a smooth curve in Q intersecting 
{Ai = 0} transversely in (0, A2), where A2 £ (A 2 , A2) and k(0) = (0, A 2 ). Assume further 

(i) U (W»(E, K (ji)) x { M }) rh U (W s (P, K (fi)) x { M }) 

(ii) W'(E,(pX)) rh W u (P,(0,\ 2 )) 

Then there is a sequence (fJ. n ), hm fj,„ = such that for all A = K(fj, n ), u> 1, there is a one-homoclinic 

n— >oo 

orbit to E. 

Assumption (i) of the lemma claims that the extended unstable and stable manifold of the equilibrium and 
of the periodic orbit, respectively, intersect transversely, while assumption (ii) claims that the stable and 
unstable manifold of the equilibrium and the periodic orbit, respectively, intersect transversely. 

Now arises the question whether all K,(p n ) lie on one continuation curve as in our motivating example - 
cf. Figure 4. In panel (i) of this figure, the black dots and squares correspond to parameter values on k 
for which a homoclinic orbit exists. The different shapes indicate that the homoclinic orbits are related to 
different EtoP cycles. Indeed, our considerations in Section 6 confirm that in the intersection of k with q, 
there exist two EtoP cycles. This feature has not been considered in panel (ii). 
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Figure 4: One-homoclinic orbits on a curve k which intersects the continuation curve of the codimension-one het- 
eroclinic orbits transversely. Panel (i) is related to the laser model (1.2). The dots and squares indicate that the 
homoclinic orbits correspond to different EtoP cycles which exist at the intersection point of n and ci,. Panel (ii) 
visualises the statement of Lemma 2.1. 



Denote the Floquet exponents of P by a s (X) < < a u (A) and assume the associated Floquet multipliers to 
be positive. Moreover, let 5 > be a sufficiently small constant and / := [— S, 6], and let S 1 := M/^ 2jr , and 
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x ~2?r V x = y mod 2ir. We introduce the so-called Fcnichcl coordinates v = (v c , v s ,v u ) G S 1 xlxl, which 
are defined in a (^-neighbourhood of the periodic orbit P. These coordinates are specially tailored to the 
hyperbolic structure of P. This is reflected by the fact that the stable manifold W S (P, A) in these coordinates 
corresponds to the set {v u = 0}, whereas {v s = 0} represents the unstable manifold W U (P, A). Further 
fixing v c = ip yields the single strong stable and strong unstable fibres W ss (P(ip, A), A) and W uu (P((p 7 A), A), 
respectively. The periodic orbit itself is given by the set {v s = 0, v u = 0}. 

Lemma 2.2. Assuming Hypotheses 2.1 (ii) and 2.2 are met, there is a smooth change of coordinates locally 
near P such that x = f(x, A) becomes 

if = l + A c (v,X)v s v u , 

v s = (a s (X) + A s (v,X))v s , (2.2) 
v u = (a a (X)+A u (v,\))v u , 

where v = (v c , v s , v u ) G S 1 x / x / and A c , A s , A u are some smooth functions in v and A with 

A l (v c , 0, 0, A) = 0, i = c,s, Ul VAeQ. 

The Fenichel coordinates were introduced at first by Fenichel in the context of slow/fast systems, [12]. In 
[13] the transformation into the Fenichel coordinates in the context of slow/fast systems is described in more 
detail. The transformation near the hyperbolic periodic orbit P, and hence the proof of Lemma 2.2, is done 
in an analogous way. The idea of the proof of Lemma 2.2 can also be found in [1]. However, note that 
in [12] and [13] the Fenichel coordinates are derived merely locally, but for our purpose we need a global 
change of coordinates with respect to the periodic orbit P. To guarantee that the Fenichel coordinates 
can be introduced along the whole periodic orbit, we exploit the assumption that both nontrivial Floquet 
multipliers are positive, Hypothesis 2.1 (ii), since this implies that the stable and unstable vector bundles of 
P arc orientablc. 




V s 

Figure 5: The cross-sections E m and E° ut . 



Next we introduce sections near P. 

S in := S 1 x {v s =8}xl, S out := S 1 x / x {/ = 5}, (2.3) 

which arc illustrated in Figure 5. Indeed, these sections are connected components of the sections E m and 
vjout f rom t ne Introduction. Further we define 

T := {(<p, A) € S 1 x Q : W S {E, A) n W uu {P{v, A), A) n E ou * / 0}. (2.4) 
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Thus r consists of all the tuples (A, <p) for which there exists a PtoE connection that contains the strong 
unstable fibre W uu (P((p, A), A) to the base point P(<p, A). 

Further, let Ur be an open neighbourhood of r in S 1 x Q. 

Hypothesis 2.3. There is a smooth function g : Ur — > I and an e > such that 

{(if, v s ,6) & W S {E, A) n : \v s \ < e, fo>, A) e U r } = {{<p, g(<p, A), S) : (y>, A) 6 t/ r }. 

As a consequence of that hypothesis, we get that T coincides with the zeros of g, cf. Figure 6: 

T:={(i P ,X)€S 1 xQ:g(ip,X)=0}. (2.5) 




Figure 6: Visualisation of Hypothesis 2.3 and its consequence: In this illustration Ai is fixed with Ai = 0. T coincides 
with the zeros of g. The graph of g is only drawn for a sample of \2-values. The dashed lines indicate a possible 
continuation of g outside ofUr- 

Figure 6 does already include some specific features of g or T, respectively, which we demand in the following 
hypothesis: 

Hypothesis 2.4. (i) There is a constant b > such that \g\ 2 (ip,X)\ > b, for all (<p,X) £ Ur- 
(ii) There is a smooth function z : S* 1 x J\ — > J2 such that T = graph z. 

As a consequence of (2.5) and Hypothesis 2.4 (ii) we find 

fl ( V ,Ai,«(vJ,Ai))=0. (2.6) 
As a transversality condition for z we assume: 
Hypothesis 2.5. z v (<p,0) = => z vv (<p, 0) ^ 0. 

Fix some A° close to zero, and let ip° be some value such that z v (Lp a , Aj) = 0. Using this, we define 
A° := z(ip°,Xl), and A := (A5,A°). Now, considering the derivatives of g(-, A?, z(-, A?)) at <p — <p° we find 
with (2.6) and Hypothesis 2.5 that 

5 ^°,A°)=0, ffw (/,A°)^0. (2.7) 

Note that graphg(-, A ) describes the stable manifold W S (E, A ) near ip°. Therefore (2.7) means that 
W S (E, A ) and W U (P, A ) have a quadratic tangency in ip°. We refer to Figure 7 for an illustration of 
the consequence of Hypothesis 2.5. 
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W U (P, A) 
r = graph z 

Figure 7: Quadratic tangency of W U (P) and W S (E) as a consequence of Hypothesis 2.5: In this illustration Ai is 
fixed with Ai = Ai. For Ai = 0, this is an enlargement of a detail around the turning point ofF in Figure 6. 



Let the function z(-,Ai) take its minimum in f m in(Xi), and similarly let z(-,Ai) be maximal in Lp max (X\). 
This defines functions Ai jTOi „ : Ai \-t z((p min (Xi), Ai) and \i, ma x ■ M ^ z((p max (Xx), Ai), both mapping 
J i —> J 2 . The graphs of these functions define the A- region for which a heteroclinic cycle connecting E and 
P exists. In our motivating example this region is just the stripe between the two curves tb - cf. Figure 2. 
Hence, the maximal interval [A 2 , A2] defined in Hypothesis 2.1 is given by 

A 2 := z(<p m i n (0),0), A 2 := z(ip max (0),0). 

Moreover, for each A between the graphs of Ai i?n m and Ai )TOax , there are at least two heteroclinic PtoE 
connection. These undergo saddle-node bifurcations on the graphs of Xi >m i n and Xi. ma x- In particular, 
moving along the A2-axis the heteroclinic PtoE connections undergo saddle- node bifurcations in A 2 and A2. 
If z has exactly one minimum (and hence one maximum), there are exactly two heteroclinic PtoE connections 
between the graphs of Xi t min and Xi^ max - cf. Figure 6. 

Next we consider the EtoP connection. 

Hypothesis 2.6. There exist smooth functions h u : Q — > I , h c : Q — > S 1 such that 

{(v c ,S,v u ) G r(E,A)nE m , A g Q} = {(v c ,S,v u ) = (h c (X)J,h u {X)), X e Q}. 

Moreover, 

(i) /i u (0, A 2 ) = 0, and VA 2 G J 2 holds h\{Q, A 2 ) ^ 0, 

(ii) 3q < 1 : G S 1 |^/i c (0,z((p,0))| < q. 



{(/ 1 c (A 1 ,A 2 ),/ l "(A 1 ,A 2 )) : Ai G 1/(0)} 




^y(hno,x 2 ),h u (o,x 2 )) W U (P) 





Figure 8: Visualisation of Hypothesis 2.6 (i) with A2 = A2: The curve (h c (\i, A2), h u (\i, A2)) = (J W U (E, Ai, A2) 
intersects W S (P) transversely. 
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By definition 



W U {E, A) n W S (P, A) n S m + <^> ft" (A) = 0. 



(2.8) 



So, Hypothesis 2.6 (i) says that for all A on the A2-axis, there is heteroclinic orbit connecting E to P. In 
other words, the A2-axis is on a par with the curve Cb of our motivating example - cf. Figure 2 or Figure 4, 
respectively Moreover, moving through the A2-axis transversely effects that the EtoP connection splits up 
with nonzero speed - cf. also Lemma 2.1. The consequences of Hypothesis 2.6 (i) for the shape and mutual 
position of the traces of W U (E) and W s (P) are depicted in Figure 8. Finally, we note that by this assumption 
h% (0, A2) is bounded away from zero. 

Recall that z(<p, 0) determines the A2 values for which a EtoP cycle exists (clearly Ai = 0), where tp is the in- 
coordinate value of the intersection of the corresponding PtoE connection with S otl *. Whereas h c (0, z(ip, 0)) 
is the ^-coordinate value of the intersection of the corresponding EtoP connection with S m . Hence, Hy- 
pothesis 2.6 (ii) yields that the proportion of the alteration rates of these -^-coordinates is bounded by q < 1. 
In other words, these coordinate values must not move against each other too fast. Despite this geometric 
interpretation, this hypothesis is more technical in nature. It will be used in the next section for solving the 
bifurcation equations. 

Now we can state our main result guaranteeing a snaking scenario. 

Theorem 2.1. Assume Hypotheses 2.1 - 2.6. Then there is a constant Lq > ; and there are functions 
Xi : (Lq,oo) — > M., i = 1,2, such that for each L > Lq there is a one-homoclinic orbit to E for A G Q that 
spends time 2L between E m and £ ou * if and only if A = (Ai(L), \2(L)). 

Moreover there are an r\ > ; a 2tt -periodic function </5g(') an d a positive bounded function a u such that 



Ax(L) = 



^ 1 (0. z (Vo(i)+2i,0)) c - 

A 2 (L) = z{&{L) + 2L, 0) + 0(e-" L ) 



(0.z( V o(L)+2L,0))L^ 1 + Qf e -nL\\ 



It follows immediately that Ai(L) tends to zero as L goes to infinity. Further it is obvious that A 2 (i) is 
a perturbation of z((Pq(L) + 2L,0). This result resembles pretty much the statement about the snaking 
parameter /1 given in [1, Theorem 2.2]. But here, in contrast to [1], the term f^L) is periodic and not 
constant. If <Pq{L) + 2L is monotonically increasing, then A2O) essentially copies the behaviour of z(-, 0), cf. 
Figure 9 and Lemma 2.3. 








A 2 


V 




h 




A(L) 




Ai 





Figure 9: The relation between graphz and the snaking curve A(L) = (Ai(L), \2(L)). The shape of graph z(-, 0) 
depicted in the left panel is passed on to the snaking curve A(L) in the right panel. The snaking curve accumulates 
at {0} x I2, the set of parameters for which the primary EtoP cycle exits, cf. Hypothesis 2.1 (in). 



In the following lemma, we describe the shape of the snaking curve (Ai(i), A2(£)) somewhat closer. We 
consider 0) as a periodic function R — > R. We denote the first and second derivative of A2 by A 2 and A' 2 ', 
respectively. 
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Lemma 2.3. Assume Hypotheses 2.1 - 2.6 with the more severe condition q < 1/2, cf. Hypothesis 2.6 (ii). 
Then $:Lh) ¥q(L) + 2L is a transformation, and for each tp with z v (tp, 0) = exists a unique L in a small 
neighbourhood of $>~ 1 (tp) such that X' 2 (L) = 0. Moreover A 2 (Z/) ^ 0. These are the only zeros of X' 2 . 

The proofs of these statements are carried out in Section 3. Prior to that, however, we give a geometric 
explanation with the help of the Figure 10. Assume that the unstable manifold of the equilibrium depends 
only on Ai, and similarly that the stable manifold of the equilibrium depends only on A 2 : W U {E, A) = 
W U {E, Ai), W S (E, A) = W S (E, A 2 ). In Figure 10 (i), we consider a fixed Poincare section of P. Fix some 
A 2 - and therefore one particular position of W s (E) - and assume that an increasing A 2 effects upward motion 
of W S (E). The bullet defines a Ai for which W S (E) and W U (E) intersect and therefore a homoclinic orbit 
to E does exist. This homoclinic orbit can be continued by moving W S (E) up and down. The corresponding 
continuation curve of the bullet in the A-space is displayed Figure 10 (ii). 




l,max 



Figure 10: The creation of a snaking curve: Panel (i) shows part of a (global) Poincare section containing both E 
and P. Panel (ii) shows the continuation curve of one-homoclinic orbits to E. 



3 Snaking analysis 

This section is devoted the proof of Theorem 2.1. A one-homoclinic orbit to E can be conceived as built of 
three pieces: an orbit segment in W U (E) running from E to S m , a solution connecting E m and E ou4 and an 
orbit segment in W S (E) running from E out to E. Let v be a solution starting in S™ and arriving after time 
2L in E out . Then v belongs to a one-homoclinic orbit to E if the following coupling conditions are fulfilled 

v(0,X) € W U (E,X) nS' m , v(2L,X) g W S (E,X) nS™'. (3.1) 

As a consequence of the following Lemma 3.1 we get that for given L there exists a unique solution v 
starting in a certain submanifold of E m and arriving after time 2L in E OMt . Afterwards, in the actual proof 
of Theorem 2.1 we use these solutions to formulate coupling equations according to (3.1). The more general 
setting of Lemma 3.1 is used in Section 4. 

Lemma 3.1 (Shilnikov problem near the periodic orbit). There is a positive constant Lq such that for all 
L > Lq, all ((f,X) G S 1 x Q, and Xs,Xu & {±1} there exists a unique solution v(t), also referred to as 
v(t,<p,X,Xs,Xu), of (2.2) with 

v s (0) = Xs S, v c (0) - tp and v u {2L) = Xu S. 
Moreover there is a positive constant r\ < min> g (g{|a ;i (A)|, a"(A)} such that 
v(0) =(tp,XsS,Xua u e- 2ali ^ L (l + 0(e-^ L ))), 

(3.2) 

v(2L) =(^ + 2L + 0(e^ L ), Xs a s e 2aS ^ L (1 + 0(e^ L )), Xu 6), 
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where a s and a u are positive functions depending on (ip, \,Xs,Xu)- Moreover, a s r u \(-, -,XsiXu) are smooth. 
For the derivatives of v holds 

D^v(2L) = (D^itp + 2L) + 0(e-^ L ), Xs D^ (a s e 2aS ^ L ) (1 + 0(e^ L )), 0). 
Here € {L, A, tp} for i = 1, . . . ,j and j e {1, 2, 3}. 

Since for fixed Xs,Xu, the functions a s and a u are denned on the compact set S 1 x Q, they are bounded 
away from zero. 

In [21], Krupa and Sandstcdc consider the Shilnikov problem in the context of slow/fast systems, where the 
slow manifold possesses a normally hyperbolic structure. Lemma 3.1 is the counterpart to [21, Theorem 4], 
and the proof proceeds in the same way as it is done there. 

Proof of Theorem 2.1 

To describe the transition from S* n to S ou * we use the function v defined by Lemma 3.1. By our choice 
of E m and T, out , cf. (2.3), this transition is determined by v(-,(p, A) := v(-,(p,X,l,l). Using the notation 
introduced in Section 2, equation (3.1) translates to 

v c (0^,X) = h c (\), (3.4) 
v u (0,<p,\) = h u (\), (3.5) 
v s (2L,tp,X) = g(v c (2L,cp,X),X). (3.6) 

Recall that we arc interested in those homoclinic orbits which are in a small neighbourhood of a heteroclinic 
cycle. These cycles are determined by h u = and 5 = 0, cf. (2.8) and (2.5). For that reason we solve (3.4) 
- (3.6) near h u = and g = 0. 

In accordance with (3.2) we find that v c (2L, <p, A) = if + 2L + 0(e~ nL ). Motivated by this equality, we 
introduce the following time transformation 

21 = 2L + 0{e^ L ). (3.7) 

Indeed, (3.7) can be solved for 

L = L*{l,X) = l + 0(e- 111 ). (3.8) 

Using this, equation (3.6) can be rewritten as v s (2l + 0(e~' nl ) 1 A) = g(ip + 21, A). Altogether, using the new 
time / and the estimates (3.2) the system (3.4) - (3.6) reads 

^ = h c (\ 1 ,\ 2 ), (3.9) 
a u e~ 2aU(x ^\l + 0(e-"')) = h u (X 1 ,X 2 ), (3.10) 
a s e 2a3 ^< x *»(l + 0(e-'")) = g(<p + 21, A x , A 2 ). (3.11) 

First, we consider (3.11), which describes the coupling in Y, out . Recall that we want to solve (3.6), and 
therefore also (3.11), near 5 = 0. Furthermore, recall from (2.6) that g(ip + 21, Ai, z((p + 21, Ai)) = 0. Now 
write 

A 2 = z{<p + 2l,\x) + fi 
and expand g(ip + 21, Ax, z(ip + 21, Ai) + /i) w.r.t. fi. Inserting in (3.11) gives 

Oae 2a*(Ai,«(*H-ai 1 A 1 )+ / 0i( 1 + 0(e-"')) = gx 2 (<P + 21, \ u z(<p + 21, Ar)) M + 0{p?). 
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Using this and Hypothesis 2.4, the coupling equation (3.6) eventually reads: 

Zt+2\t^^ + = f + <V)- ( 3 - 12 ) 

For <C 1. sufficiently large I and all ip this equation can be solved for p, = p?(l, <p, Ai) by means of the 
implicit function theorem [8, Chapter 2.2]. The solving function ji* is differentiable. Further we see from 
(3.12) that 

Altogether we find that the coupling equation (3.6) can be solved for A2 = A2^, <p, Ax) with 

X 2 {l,p, Ai) = z(ip + 21, Ai)) + n*(l,(p, Ai), 
where the leading order term of p* is given by (3.13). 

Now we turn towards the coupling in £ m which is determined by (3.9) and (3.10). Using the representation 
of A2 , these equations read 

if = h c (Xi, z{<p + 21, Ax) + p*(l, if, Xi)), (3.14) 
aue -*<* u {^A<p+UM)+f{i,vM))l(i + 0(e~ 7 ' 1 )) = h u (Xi,z(ip + 21, Ai) + p*{l, ip, Ai)). (3.15) 

We solve (3.14), (3.15) for (<p, Ai) depending on I. Note that <p is the v c coordinate where the prospective 
homoclinic orbits hits E m . Hence, ip may vary within a "large" range. To handle this difficulty analytically, 
we first consider the "unperturbed equation" tp = h c (0,z(ip + 21,0)). For that we consider z(-,0) as a 27T- 
periodic function R — > R, cf. Hypotheses 2.1 and 2.4. Hence, h c (0,z(- + 21,0)) is a 27r-periodic function 
as well. Because of Hypothesis 2.6 (ii), we can apply again the implicit function theorem to find a unique 
solution (Pq(1) on R such that 

p = h c {0, z{tp + 21, 0)) (p = y* {l). (3.16) 
Note that <Pg(-) is again 27r-periodic. Now, write 

<p = ip* Q (l) + ip, 

and we define H c (l, ip, Ax), H u (l, tp, Ai) by 

H°/ U (l, V, Ai) := h c ' u (\uz(<pZ(l) +4> + 2l, Ax)). 
Using these terms the right-hand sides of (3.14) and (3.15) read 

h c '»(\ lt z(<p* (!) +1/J + 21, Ax) + n* {I, <p*(l) + i>, Ax)) = H c ' u (l, i/>, Ax) + W>*(Z, <p' (l) + i/>, Ax)), 
where r c / u (p*) = 0(p*). Further, since ip* is the unique solution of (3.16), we find 

H c (l,0,0) = h c {0,z(<p*(l) + 21,0)) = <p* {l). 
In accordance with Hypothesis 2.6 (i), we find furthermore 

H u (i,o,o) = o, h;(i,o,o) = o. 

Hence, (3.14) and (3.15) are equivalent to 

rp = H$(l, 0, 0)i/j + H c Xl (1, 0, 0)Ax + 0(|(V, Ax)| 2 ) + r c (p*), (3.17) 
„ e -!»-(W,AOi( 1 + ( e -fl«)) =fl 5. i (| ) o j o)Ai + 0(|(^Ai)| a )+r w O**) ) (3.18) 

with a u {l, ip, Ax) := a"(Ax, z(p* (l) + ip + 21, A x ) + p*(l, <p* {l) + i/>, Ax)). 
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Our goal is now to solve the system (3.17), (3.18) for (i/>,Ai) depending on I. To this end we invoke again 
the implicit function theorem. The main observation is that due to Hypothesis 2.6 

|J35(i,0,0)|< 9 <l, #^,0,0)^0. 

Note that due to (3.3), the corresponding partial derivatives of r c (fj,*(l, <Pq(1) + tp, Ai)) and r u (fi*(l, <fo(l) + 
tp,Xi)) tend to zero as I — > oo. Therefore there exist unique functions ip*(l), A*(Z) satisfying the system 
(3.17), (3.18). Accordingly (3.14), (3.15) are satisfied by 

0*(O:=Vo(O + ^*(O and A i(0- 
Inspecting (3.15) and (3.17), we find with a u (l) := a u ((po(l), 0, z((Pq(1), 0)) that 

A *(Z) = g „ a " ( g, 0) e- 2a "('- -°)'(l + 0(e-" ! )) and = 0(e~ 2&u ^ 1 ). 

Altogether, for (3.9) - (3.11) we find the unique solution (tp, Ai,A 2 )(Z) = (<t>*(l), X*(l), X%(1)), where 

A^(/) := A 2 (Z, 0*(Z), AKO) = *#*(0 + 2/, A*(0) + <f (0, A*(Q) 

= z(^(Z) + 2Z,0) + O(e-" i ). 

Note that v spends time 21 + 0(e~ vl ) between E m and T, out , cf. (3.7). So, in view of the statement in 
Theorem 2.1, we define 

4>{L):=4>*{l{L)), Ai(L) := A 1 (Z(i)), A 2 (L) := AJ(i(L)). 
Then (<£, Ai, A 2 )(Z) = {<j>{L), Xi(L), A 2 (L)) solves (3.4) - (3.6). The above considerations yield 

0(i)=^(L) + O( e -' L ), 

AiW - 7?0^)) e " 2<i " (i ' O ' O)L ( 1 + O^))' 
A 2 (L) = z(^(L) + 2L, 0) + 0(e-" L ). 
This finally completes the proof of Theorem 2.1. ■ 

Proof of Lemma 2.3 

In what follows, we sketch the proof of Lemma 2.3. We note that, due to (3.3), the O-term in the repre- 
sentation of A 2 (L) is differcntiablc and its derivative can be estimated by a O-term of the same order. The 
same holds true for higher derivatives. Therefore we find 

X' 2 (L) = z v (tp* (L) + 2L, 0)«(L) + 2) + 0(e^ L ). 

Recall the determining equation <fio(L) = h c (0, z(<Pq(L) + 21,0)) for ipq(L), cf. (3.16). From this equation, 
we get an estimate of the derivative of <Pq(L), whereby we finally confirm that for sufficiently large L 

V l\L)+2^0. 

So, necessarily the zeros of X' 2 (L) are close to the zeros of z v ((po(L) +2L, 0)). Let ^(^(Lo) + 2L ,0)) = 0. 
Using the contraction principle we find a neighbourhood U(Lq) of L$ in which A' 2 (L) = has a unique 
solution L. Straightforward computations show A 2 '(L) ^ 0. Further, the size of the neighbourhood U(Lq) 
can be chosen independently of Lq. Outside the union of these neighbourhoods, z v {-, 0) is bounded away from 
zero. This finally shows that outside the union of these neighbourhoods A 2 (L) has no zeros for sufficiently 
large L. ■ 
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4 Negative Floquet multipliers 



Now we discuss the scenario, where the nontrivial Floquet multipliers of the periodic orbit P are negative - 
in other words, we replace Hypothesis 2.1 by: 

Hypothesis 4.1. The nontrivial Floquet multipliers of P are negative. 

Recall, since our setting is in R 3 , the two nontrivial Floquet multipliers of P must have the same sign. 
Negative multipliers cause that the vector bundle consisting of the eigenvectors of the monodromy matrices 
along P is a Mobius strip and thus not orientable. Hence, we cannot introduce Fenichel coordinates near P. 
We overcome this difficulty by introducing local coordinates, which are not 2tt periodic, but of period 47r. 
To this end we transform at first (2.1) into normal form, cf. [26, Theorem 3.11], which gives: 

6 = 1, 

y = B(6)y + F(6,y,\), 

where y = (2/1,2/2) £ K 2 and 9 G S 1 ; B and F are smooth. Furthermore we straighten the stable and unstable 
fibres of P, as it is done in [13]. This yields that the function F satisfies 

F{6, 0, 2/2, A) = and F(8, 2/1, 0, A) = and D y F(9, 0, 0, A) = 0, VAgQ, V6> G S 1 . 

After that we apply Floquet theory to the linear system y = B(9)y, see [26, Theorem 3.12] for more details. 
This transforms the above normal form into: 

v c = 1, 

V s = a s (X)v s +F s (v c ,v s ,v u ,X), (4.1) 
v u = a u (X)v u +F u (v c ,v s ,v u ,X), 

where v c G S 1 := R/^ 4?r , and x ~4 W y <f> x = y mod47r. Moreover, similar arguments as leading to 
Lemma 2.2 yield F s (v c , 0, v u , A) = F u (v c ,v s ,0, X) = 0, as well as, D v sF s (v c , 0, 0, A) = D v uF u (v°, 0, 0, A) = 0. 

Note that by this construction the two points (v c , v s ,v u ) and (v c + 2tt, —v s , —v u ) represent the same point 
in {9, 2/)-coordinates. In other words, two points are identified via the map: 

i: S 1 x I x I ^S 1 x I x I 

(v c , v s , v u ) ^ (v c + 2tt, -v s , -v u ). 

Next, in accordance with the procedure in Section 2, we introduce a cross-section S m of W S (P) intersecting 
orthogonally the stable fibres of P in a distance 6 of P. Similarly we define E out . In (v c , v 3 , w")-coordinates 
these sections read: 

:= S 1 x {v s =6}xl, Y,° + Ut := S 1 x I x {v u = 6}. (4.2) 

The subscript "+" refers to the positive value S for the fixed v s - and ^"-coordinate, respectively. The sections 
defined in (4.2) are identified via the map i with 

S m := S 1 x {v s = -5} x /, £ rat := S 1 x I x {v u = -5}. 

Further we introduce a set T similarly to (2.4) - formally replacing 5* 1 by S 1 : 

r := {{if, A) G S 1 x Q : W S (E, A) n W uu {P{<p, A), A) n Z° + Ut / 0}. (4.3) 

Note that each strong unstable fibre W uu (P(ip, A), A) of P intersects E out twice. In the terminology of (4.3) 
those two points arc represented by W uu (P((p, A), A) n S^ Mt and W uu (P(ip + 2ir, A), A) n S^ 1 '. 

Let Ur be an open neighbourhood of T in S 1 x Q. Regarding the PtoE connecting orbit we assume: 
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Hypothesis 4.2. There is a smooth function g : Ur — > / and an e > such that 

{(if, v s , 6) G A) n E?* : H < e, {<p, A) G C/ r } = g&, A), 5) : {<p, A) 6 C/ r }. 

Hypothesis 4.3. (i) There is a constant b > smc/i that \g\ 2 (<p, A) > 6, /or aZZ (93, A) G Ur- 
(ii) There is a smooth function z : 8 1 x J\ — > J 2 such that T = graph z. 

Consequently 

g{vAi,z{ip,\i)) = 0. 

Hypothesis 4.4. z v (ip,Q) = z vv (<p,0) ^ 0. 

Finally, regarding the EtoP connecting orbit we assume: 

Hypothesis 4.5. There exist smooth functions h u : Q — > I, h c : Q — > S 1 such that 

{(v c ,S,v u ) G VF U (£,A) nS™ A G Q} = {(w c ,<5,v u ) = (h c (X),S,h u (X)), A G Q}. 

Moreover, 

(i) h u (0, A 2 ) = 0, and VA 2 G J 2 holds h^(0, A 2 ) > 0, 
(ii) 3q< 1 : V^GS 1 1^(0,^,0))! <<?. 

Indeed, in Hypothesis 4.5(i) it already suffices to assume /i" (0, A2) 7^ 0. The specification stated in the 
hypothesis determines the sign of the functions A^~ and Aj~ in the way as stated in the theorem below. 

Now, the analogue of Theorem 2.1 reads: 

Theorem 4.1. Assume Hypothesis 2.1, and Hypotheses 4-1 - 4-5. Then there is a constant Lq > ; and 
there are functions Xf , A~ : (Lq, 00) — > R, i = 1, 2, such that for each L > Lq there is a one-homoclinic orbit 
to E for X G Q that spends time 2L between E m and £ out if and only if A = (Xf'(L), Xf(L)). 

Moreover there are an r\ > 0, two Att -periodic functions <Pq(-), V?o~(') an d positive bounded functions a+, a~ 
such that 

xf{L) = 

Xf(L) = z(<pft(L) + 2L, 0) + 0(e-" L ). 



A visualisation of the statement of this theorem is given in Figure 11. 





A-(L) 




A+(L) 



Ai 



Figure 11: As in the case of positive multipliers the shape of graph z(-, 0) is passed on to the snaking curves A ± (L), 
cf. Figure 9. The snaking curves accumulate at {0} x I2 from different sides. 
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Proof. We pursue the same strategy as in the proof of Theorem 2.1: We construct one-homoclinic orbits 
to E by coupling in E m and E out the unstable and stable manifolds, respectively, with solutions according 
to Lemma 3.1. However, here there are two possibilities for the transition from E m to E otlt . There are 
solutions of (4.1) starting in E™ (~l {v u > 0} and end up in ET"' n {v s > 0}, and there are solutions starting 
in E™ fl {v u > 0} and end up in E°_ u * n {v s < 0}, cf. Figure 12. In the language of Lemma 3.1 this distinction 
is determined by the signs of Xs and Xu- 




Hypothesis 4.5. In E+ 1 there is drawn the graph of g(-,X) for fixed X, cf. Hypothesis4-2. 



The transition E™ n {v u > 0} to E° ut n {v s > 0} 

Here we employ solutions of Lemma 3.1 with (x s , Xu) = (1, 1) what we suppress from the notation. In this 
case the argumentation runs completely parallel to the proof of Theorem 2.1. We confine to sketch the 
procedure: The coupling equations analogue to (3.4)- (3.6) are almost the same: 

v u (0,<p,\)=h u (\), 
« S (2L,^,A)= 5 K(2L^,A),A), 

with the only difference that here ip,v c 6 S 1 . From that we gain the analogue to (3.9)-(3.11) 

V = h c (\), 

a u e- 2aU{x)l {\ + 0(e-"')) = h u (X), (4.4) 
a s e 2a 'W l (l + 0{e- r ' 1 )) = g{<p + 21, A). 

Proceeding in the same way as in the proof of Theorem 2.1 we get the solutions \ f(L) and A^(L) as stated 
in the theorem. The function tp^ solves the analogue to (3.16). 

The transition E m n {v u > 0} to n {v s < 0} 

First we note that the intersections of W U (E, A) with E m written in the form (h c (X),S,h u (X)) e E™, cf. 
Hypothesis 4.5, arc identified with (h c (X) + 2ir, — 6, — h u (X)) G E" 1 . This allows to employ solutions of 
Lemma 3.1 with (xstXu) = ( — 1, 1) for our analysis. Again we suppress the Xs , x« _ dependence from the 
notation. 



18 



Thus the coupling equations read: 



v c (0 7 tp,\) = h c {X)+2n 7 
v u (0,<p,X) = -h u (X), 
v s (2L,<p,\)=g(v c (2L,tp,\),V. 

With the results of Lemma 3.1 this can be rewritten as 

tp = h c (\)+2n, 



ue -aa«(A)i( 1 + {e~ r ' 1 )) = -h u (X), (4.5) 
2 « a W(l + 0(e-"')) = g(f + 2Z, A). 



a 

-a.,e 



Now we can proceed again as in the proof of Theorem 2.1, with the minor difference that the function ip 
results from the fixed point equation 

tp = h c (0,z(tp + 2l,0)) + 2n. 
This finally leads to the solutions X^(L) and A^"(L). 

We want to note that A^ and A^f have different signs. This follows immediately from (4.4) and (4.5) in 
combination with Hypothesis 4.5(i). I 



5 A nonsnaking scenario 

Now we consider a further continuation scenario of homoclinic orbits. We retain the general Hypothesis 2.1, 
and assume positive Floquet multipliers. Further we adopt the definition of T, cf. (2.4), and assume 
Hypothesis 2.3. In Section 3, we have seen that the shape of the continuation curve (Ai(-L), A2(£)) is basically 
determined by the form of the set T, which was assumed to be the graph of a function z : S 1 x Ji — >• J 2 . Now 
we investigate the consequences of altering the corresponding Hypothesis 2.4. More precisely, we suppose 
that r is no longer graph of a function, but for fixed Ai a closed curve in each case. 

In contrast to the situation in Section 3, we end up with a sequence of closed continuation curves, so-called 
isolas. Consequently, the one-homoclinic orbits close to the primary EtoP cycle do not lie on one (global) 
continuation curve. The addressed curves tend to {0} x I2 in the sense of the Hausdorff metric . 

A similar scenario was already discussed in [1] for the Hamiltonian case under slightly different assumptions. 

We assume the following: 

Hypothesis 5.1. Let I\ C J\ be a closed interval containing zero and let I v C S . There exist smooth 
functions tp : S 1 x 1\ — > I v and A2 : S 1 x I\ — > Ji such that 

T = {(0(r,X 1 ), Ai, A 2 (r,Ai)) : r e S 1 , Aj e h}, 

where 

(D r <p(r,\ 1 ),D r X 2 (r,X 1 ))^0, Vr g S 1 , V Ai 6 h- 

Hence, the set Tn {Ai = const.} is a closed, regular curve in I v x J2 parametrised by some parameter r € S 1 . 

Again we denote by Ur an open neighbourhood of r in S 1 x Q. Further, (2.5) together with Hypothesis 5.1 
yield as counterpart to (2.6) 

<?(0(r,Ai),Ai,A 2 (r,Ai)) = O, V(r, Ai) G S 1 x h. (5.1) 
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Figure 13: Visualisation of Hypothesis 5.1. Depicted is the closed curve {(if>(r,0), Ai, \2(r,0)) : r G S 1 } = {<?(■, 0, •) = 
0}. 

Hypothesis 5.2. 

D v g($, 0, X 2 )D r X 2 {r, 0) - Dx 2 g(0, 0, X 2 )D r $ 2 (r, 0) + 0. 

Assume Hypothesis 5.1. Then Hypothesis 5.2 means that the gradient of g(-,0, •) docs not vanish at any 
point within the set T. 

Hypothesis 5.3. 

{<p(r,Q) ~ h c (0,X 2 (r,0)), r e S 1 } C S 1 . 

The subtraction in the hypothesis is done in S 1 . Since l v is a proper subset of S 1 , this hypothesis is satisfied 
if | J D A2 ^ c (0,A 2 (r,0))^ r A2(r-,0)| «l,Vr£ S 1 . 

Theorem 5.1. Assume Hypotheses 2.1-2.3, 2.6 (i) and Hypotheses 5.1-5.3. Then, there is a sequence of 
mutually disjoint closed continuation curves Afc ■= {(Ai.fc(r), X 2 .k(r)), r & S 1 }, k G N. These curves tend 
towards a segment of the X 2 -axis in the sense of the Hausdorff metric. More precisely, for sufficiently large 
k G N, there exist mutually disjoint intervals Ik C R and smooth functions Xi.k : S 1 — > Ji, i = 1,2, and 
Lk ■ S 1 — > Ik such that there is a homoclinic orbit to E for A G Q with flight time L G Ik from E m to T, out , 
if and only if, there exists an r G S 1 such that A = (Ai^r), A2 j / C (r)) and L = Lk(r). Moreover 

A life (r)=0( e -" fc ) 

A 2 , fc (r) = X 2 (r,X h k(r)) + 0(e^ k ). 



Afc A fe+1 

OOOl 


A 2 
h 


Ai 





Figure 14: Visualisation of Theorem 5.1: The continuation curves Afc. 
Proof of Theorem 5.1 

We follow the lines of the proof of Theorem 2.1 up to equations (3.9) - (3.11), which we repeat here: 

<p = h c (X 1 ,X 2 ), (5.2) 
au e-2«"(AiA.)i( 1 + ( e -tf)) =h u (X 1 ,X 2 ), (5.3) 
a s e 2a "^^ l (l + 0(e-"<)) = gfa + 21, Ai, A 2 ). (5.4) 
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Again wc look for solutions of these equations near the set T. For sufficiently small Ai £ J\ each point of a 
small tubular neighbourhood of L n {Ai = Ai} has a unique representation: 

(ip + 21) mod 2ir = ip(r, Ai) + (p 

A 2 = A 2 (r, Ai) + A 2 , 

where (<p, A 2 ) is in the normal space {(D r ip(r, Ai), D r X 2 (r, Ai))}- 1 of T at the point (<p(r, Ai), A 2 (r, Ai)). 
Hence, <p and A 2 have to satisfy the additional equation: 

D r ip(r, \ ± ) ip + D r \ 2 {r, \ ± ) \ 2 = 0. (5.5) 

At first we solve equation (5.4) together with equation (5.5). Expanding g w.r.t. (ip, A 2 ) we obtain 

g{tp + 21, Ai, A 2 ) = g(0, Ai, A 2 ) + D^p, Ai, A 2 )<£ + D X2 g(X 2 , Ai, A 2 )A 2 + 0(|(& A 2 )| 2 ). 

Inserting this into (5.4) (and in the process exploiting (5.1)) yields 

£>„ff(& Ai, A 2 )^ + ^(0, Ai, A 2 )A 2 + 0(|(fc A 2 )| 2 ) - ^'(^(nAO+A,)^ + (e^ 1 )) = 0. (5.6) 

Because of Hypothesis 5.2, we may apply the implicit function theorem to the system ((5.5), (5.6)), and find 
a unique solution ((p* , A|)(i, r, Ai) for sufficiently small Ai and \{(p, A 2 )| and sufficiently large I. This solution 
tends to zero uniformly in r and Ai as I tends to infinity. Moreover, {<f>* , A^)(^, r, \\) is diffcrentiable w.r.t. 
I and Di(tp*, X^il, r, Ai) = 0(e~ vl ). The latter can be seen by differentiating (5.6) and taking Lemma 3.1 
into consideration. 

Next we consider (5.3), where we insert ((/?*, A£). Due to Hypothesis 2.6 (i), the resulting equation can be 
solved in the same way as (3.10) in Section 3. We obtain the solution Ai = Ai(Z, r), and we find that both 
Xi(l,r) and Di\i(l,r) are of order 0(e~ vl ). The estimates for the derivatives follow from Lemma 3.1. 

It remains to solve equation (5.2), which can be written as 

(<p(r, Ai(i, r)) + 0(1, r) - 21) mod 2tt = h c (X l (I, r), X 2 (r, Aj (I, r)) + X 2 (l, r)), (5.7) 

where (<p, X 2 )(l, r) := ((p* , X 2 )(l, r, A*). To solve (5.7), we have to overcome similar obstacles as when solving 
equation (3.14). For that purpose we define 

2l (r) :=0(r,O)-h c (O,X 2 (r,O)). (5.8) 

Note that l (r) : S 1 -> R is smooth. Now we set 21 = 2l (r) +21 + 2kn, for some k £ N. Fixing k € N, we 
define 

Ai, fe ([,r) := Ai(/ (r) +Z~+fc7r,r), A 2 . fc (f,r) := A 2 (Z (r) + 1 + kir, r), <p k (l r ) ■= <PQo(r) + I + kn,r). 
Using this, we rewrite (5.7) as 

0fc(r,r) + O(Ai,fc) + O(Aa lfc ) = 21 (5.9) 

This equation can be solved by means of the contraction principle. For that we note that the terms on the 
left-hand side together with their derivatives are of order 0(e~ vk ). Thus, for each fixed r and sufficiently 
large k, equation (5.9) possesses a unique fixed point ll(r). Moreover, 1% : S 1 — ^ R is smooth, and l%(r) is of 
order 0(e~ r,k ). All in all we obtain the unique solutions 

h(r) = lo(r)+h(r)+k7T, (5.10) 
Ai, fc (r) = Ai(i fc (r),r), 
Aa,fc(r) = A 2 (r,Ai I fc(r)) + A2G i b(r),r). 
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Obviously, A& := {(Ai^r), A 2 ,fc(r)),r G S 1 } arc closed curves. Finally, with the transformation 

L fc (r) = Z fc (r) + 0(e-*«), (5.11) 



cf. (3.8), we get 

Xi, k (r) = MLk{r) + 0(e- r > lk ^),r) = 0( e -" fc ), 

A a ,*(r) = A 2 (r, A 1)fe (r)) + A 2 (L fe (r) + 0( e -*M),r) = A 2 (r, Ai, fe (r)) + 0( e -" fc ). 

Furthermore, we define the intervals = [/ fe ,7fc] := L^iS 1 ). Due to Hypothesis 5.3 and (5.8), the length of 
the interval loiS 1 ) is less than ir. Finally, since h(r) = 0{e-^ k ), it follows with (5.10) and (5.11) that there 
is a d > such that for sufficiently large k 

l_ k+1 -h>d. (5.12) 

Hence, for sufficiently large k and k 

T k C\l' k = 0. 

It remains to show that the curves h k are mutually disjoint: Assume that there exist k, k G N, k > k and 
r,f G S 1 such that 

(Ai,A 2 ) := (A 1)fc (r),A 2)fc (r)) = (A 1>fc (f), A 2 ^(f)). 
Hence, /i"(Ai, A 2 ) = h u (Xi.k{r), A 2l fe(r)) = /i"(A 1 g(f), A 2 and from (5.3) we deduce 

e -2a"(A 1 ,A 2 )( fc ( 1 + ( e -*)) = g-aa-tAi.^,)!^! + 0(6"^)), (5.13) 

where we exploited the fact that a u depends only on Ai, A 2 and ip, together with equation (5.2). Because of 
h( f ) - l k( r ) > d > cf - ( 5 - 12 ). wc infcr from ( 5 - 13 ) that 

(1 + 0(e-*)) = e -^(\iM)[l k 'l k \(i + 0(e-^)) < e -2«"(Ai,A 2 )d + e -2a«(X 1 ,A 2 )d ( e -^^_ (5.14) 

Taking the limit k, k — > oo, we see that the left-hand side of (5.14) tends to 1. On the other hand the 
right-hand side of (5.14) is close to e ~ 2aU ( x iM)<i < \ This yields a contradiction, and for this reason the 
curves Afc and Ar cannot intersect. ■ 



6 Numerical verification of the hypotheses 

In this section we show numerically that the assumptions stated in section 2 are satisfied for our motivating 
example (1.2). Recall that the setting of system (1.2) is exactly opposite to the setting described in the 
hypotheses in Section 2, i.e. dimW s (E) = 1, meaning that the orbit 7 PtoE is of codimension one and the 
orbit 7 EtoP is of codimension zero. In order to adapt the following numerical computations to the explanations 
given in the preceding sections we consider the vector field —F(x, y, tp, v\, z/ 2 ). Further wc define 

u:=(x,y,tp), A:=(^i,i/ 2 ), f(u, A) := —F(u, A). 

The statements of Hypothesis 2.6 are trivially satisfied: As shown in [20], the EtoP connection which exists 
for parameter values on cf. Figure 2, is detected by constructing a numerical test function that is defined 
by a signed distance of the endpoint of an orbit segment and the starting point of another orbit segment. 
Each zero of the test function then corresponds to a real EtoP connecting orbit. It turns out that the roots 
of the test function are indeed regular, which means that Hypothesis 2.6 is numerically satisfied. What 
is more, the way the computations are set up immediately provide the Floquet multipliers of the periodic 
orbit, hence we can easily check that Hypothesis 2.1 is also satisfied (see also below for the computation of 
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the Floquct multipliers). In the following we focus on the numerical verification of Hypotheses 2.4 and 2.5, 
which are not trivial to check. 

In our further considerations we restrict ourselves to verify the shape of T rather than to present a numerical 
verification of all requirements stated in Hypotheses 2.3 - 2.5. However, the computation of the shape of 
r requires advanced numerical techniques. In the following section, we show how to do this by finding 
and continuing the codimension-0 PtoE heteroclinic connections 7 PtoE for parameter values (^1,^2) along 
the curve Cf, (between the two intersection points with to, cf. Figure 2). More precisely, we compute the 
intersection points of 7 PtoE with a torus T — P x S* 1 surrounding P and transform the coordinates of these 
intersection points to the required format. Note that this means that T appears as a curve similar to the 
sketch in Figure 6. 

The implementation of the method 

In the following, we explain the different continuation runs (i.e. numerical solutions of boundary value prob- 
lems with varying parameters) needed for the computation of the curve T as defined in (2.4) for system (1.2). 
For the actual computations, we utilise the software package AUTO [10], which requires us to formulate the 
vector fields in the time-rescaled form u = Tf(u, A), T £ [0, 1]. 

Step 1: Finding 7 PtoE 

Similar to the computation in [20], it is possible to find the codimcnsion-0 PtoE connecting orbit 7 PtoE by 
setting up and solving appropriate boundary value problems. 

We start by constructing orbit segments u + C W S (E) and u~ C W U (P) (both reasonably close to 7 PtoE ) 
such that u + (0),u~(l) £ E, where E is a cross-section of 7 PtoE dividing the phase-space such that E and P 
are separated. Then we close the gap u + (0) — u~ (1), which corresponds to finding a numerical representation 
of 7 PtoE - Note that this methods allows us to find both possible PtoE connections. 

In the following, the boundary value problems that are used for the consecutive continuation runs with 
AUTO are listed. Note that the equilibrium point E as well as its eigenspaces are analytically known and 
hence are omitted in this listing. 

In order to compute the orbit segments, we need to continue several objects simultaneously. This can be 
achieved by extending the system by the additional objects in consecutive continuation runs. 

We assume that we have prepared the system such that A = (^1,^2) £ Cf,, half way between the two 
intersection points with Moreover, we assume that the numerical representation of P and the numerical 
representation of the unstable Floquct bundle along P and the value of the unstable Floquet multiplier, 
which we here denote by \i, is known. For more information about the necessary steps to achieve this, we 
refer to [20]. 

The first object we need to continue is the numerical representation up of the periodic orbit P. The following 
standard boundary value problem (BVP) for periodic orbits is used: 



Jo 

Note that during the continuation run, up denotes a solution of up from the previous computation step. 



u p 







Tf(u P ,X), 
up(0)-u P (l), 



(6.1) 
(6.2) 




(6.3) 
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In addition to up, the Floquct bundle up is continued using: 

u F = TD u f{u P (t), X)u F {t) + In (6.4) 
= sgn(/i)u F (0)-w F (l), (6.5) 
= {u F (0),u F (0)) -1- (6.6) 

Note that we assume that the Floquct bundle is normalised such that 1 1 up (0)|| = 1. For more information 
on this BVP, sec [11]. 

The orbit segment u~ C W U (P) is continued using the following equations: 

«-=T-/(«- A), (6.7) 
= u P (0) + Su F (0) -u-(0), (6.8) 
0=(u _ (l)-o-,n s )-J7 _ . (6.9) 

Here, a is an arbitrary point in E and ris is a normal of the cross-section E. The parameter S is initialised 
with — 10~ 4 , the solution u~ is initialised with a constant value of itp(0) + <5up(0). 

The first continuation run uses the ODEs (6.1), (6.4) and (6.7). The boundary conditions are (6.2), (6.5), 
(6.6), (6.8), (6.9), the integral condition is (6.3). The primary continuation parameter is (initialised with 
0), the remaining continuation parameters are /x, T and rf~ . The value of rj~ is considered a signed distance 
of the endpoint u~(l) to E, hence if rf~ = the orbit segment u~ ends in E and the continuation run stops. 

The second orbit segment u + C W S (E) is defined by the following equations: 

u+ = T+f(u+,\), (6.10) 
= E + ev s - u + (l), (6.11) 
0=<u + (0)-o-,n s )-7? + . (6.12) 

Here, v s denotes a vector within the two-dimensional stable eigenspace at E, which is analytically known. 
The parameter e is initialised with 10 -6 , the solution u + is initialised with a constant value of E + ev s . Note 
that we keep the boundary condition (6.11) throughout all following continuation runs, which is only valid 
because E is a saddle-focus. 

The second continuation run uses the ODEs (6.1), (6.4), (6.7) and (6.10). The boundary conditions are (6.2), 
(6.5), (6.6), (6.8), (6.9), (6.11), (6.12), the integral condition is (6.3). The primary continuation parameter 
is T + (initialised with 0), the remaining continuation parameters are /z, T, T~ and ij + . The value of r/ + is 
considered as a signed distance of the starting point u + (0) to E, hence if i] + = the orbit segment u + starts 
in E and the continuation run stops. 

After these two continuation runs, the orbit segment u~ ends in E and u + starts in E. We define z := 
(w + (0) — it~(l))/||u + (0) — u~(l)|| and initialise the new parameter r\ with the value ||u + (0) — 

Then we replace the boundary condition (6.12) by 

= u + (0)-u-(l)-r}z, (6.13) 

which means that we force the difference it + (0) — u~(l) to be in the linear subspace defined by the vector z, 
while the parameter 77 measures the gap between these two points. 

In order to close the gap, we perform a third continuation run using the ODEs (6.1), (6.4), (6.7) and (6.10). 
The boundary conditions are (6.2), (6.5), (6.6), (6.8), (6.9), (6.11), (6.13), the integral condition is (6.3). 
The primary continuation parameter is 77, the remaining continuation parameters are fi, T, T~ and T + . The 
value of 77 is considered as a signed distance between the two points u + (0) and u _ (l), hence finding a root of 
77 corresponds to finding a numerical representation of 7 PtoE . Note that there are two incarnations of 7 PtoE , 
which correspond to two different roots of 77. This concludes the first step of our method. 
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Step 2: Computing the intersection point 7 PtoE n T 

In this step we compute the intersection point of 7 PtoE with the torus T. We use the resulting orbit segments 
from step 1, but we no longer force the points u~(l) and u + (0) to be in S, instead we let u~(l) vary along 

7PtoB- 

For technical reasons, we need to include a second copy of P into the continuation using the equations: 

h P =ff(u P ,X), (6.14) 
= u P (0) -up(l), (6.15) 

= / (il P (T),U P {T))dT. (6.16) 

Jo 

The solution u is initialised with the solution u Pl and the parameter T is initialised with the parameter T. 
Recall that the parameters T, T, T~ and T + are necessary to transform the vector field such that the time 
interval on which the orbits are computed is always [0,1]. 

The following boundary condition is used to measure the distance of u~(l) to the Poincare section of P at 
fip(0): 

= (u P (0) - «-(i), /(Mo), A)) - m- (6-17) 

The fourth continuation run uses the ODEs (6.1), (6.4), (6.7), (6.10) and (6.14). The boundary conditions 
are (6.2), (6.5), (6.6), (6.8), (6.11), (6.13), (6.15), (6.17), the integral conditions are (6.3) and (6.16). The 
primary continuation parameter is rji, the remaining continuation parameters are pi, T, T, T + , T~, 5, e. 
Any root of r\\ means that the point u~(l) = u + (0) is in the Poincare section defined by up(0) and having 
the normal /(wp(0),A). Note that in general there may be several roots of rji. We choose the solution for 
which u~(l) is in the unstable fibre of up(0), i.e. u~(l) € W uu (up(0), A). In practise, we achieve this by 
using the solution for which the distance ||u _ (l) — up(0)|| is smallest. 

In the final continuation run for the detection of the intersection point, we measure the distance of u~ (I) to 
up(0) using the boundary condition 

0=||fip(0)-u"(l)l|-»72- (6.18) 

The fifth continuation run uses the ODEs (6.1), (6.4), (6.7), (6.10) and (6.14). The boundary conditions 
arc (6.2), (6.5), (6.6), (6.8), (6.11), (6.13), (6.15), (6.17), (6.18), the integral condition is (6.3). Note that 
we omit the integral phase condition of up. The primary continuation parameter is 772, the remaining 
continuation parameters are fi, T, T, T + , T~, 5, e. When the value of 772 reaches the desired radius of T, 
the continuation run stops and the intersection point is found. The actual computations are performed for 
the radius 0.1. 

Step 3: Continuation of 7 PtoE along q, 

Instead of using the (so far fixed) parameters v\ and v-i directly, we define a smooth function A : [0, 1] — > IR 2 , 
s n- A(s), such that A(0) = c& PI tb (the lower intersection) and A(l) = q, n tb (the upper intersection) 
and A([0, 1]) C Cb- Using this definition of A, we can continue the system consisting of (6.1), (6.4), (6.7), 
(6.10) and (6.14), with boundary conditions (6.2), (6.5), (6.6), (6.8), (6.11), (6.13), (6.15), (6.17), (6.18) and 
with integral condition (6.3). The primary continuation parameter is s (initialised to 0.5), the remaining 
continuation parameters are /j, T, T, T + , T~, 5, e. This final continuation run is performed for increasing 
s until a limit point for s = 1 is reached, and then for decreasing s until a limit point for s = is reached. 
These two limit points correspond to the two intersection points of Cf, with the two branches of tb. 
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7C/2 71 37C/2 (f 2K 

Figure 15: The computed curve T for system (1.2). The curve is shown as an ip-vs-V2 plot, where the angle <p is 
calculated as the local coordinate on the torus surface T and corresponds to the angle tp used in Hypothesis 2. 4- The 
curve F is interpolated between the computed points (the marked points onT) using cubic splines and it is shifted such 
that the minimum value is at ip = 0. 

This concludes our method of computing V for system (1.2), the result after transforming the coordinates of 
u~(l) = u + (Q) to the local coordinates of the torus is shown in Figure 15. Note that the shape of T verifies 
2.4(h) and 2.5, but it does not yet cover 2.3(i). In order to verify 2.4(i), the curves corresponding to graph g 
for different values of ip (cf. Figure 8) needs to be computed. The computations are very similar to the 
computations presented above and are therefore omitted here. 

7 Discussion and conclusion 

In this paper we have considered the behaviour of one-homoclinic orbits near EtoP cycles. We have rigorously 
shown how the behaviour of the primary cycle determines the continuation behaviour of the homoclinic orbits. 
Our analysis is restricted to M 3 . In higher dimcnsionsal state space the dimensions of £ m and S°™' increase 
accordingly and therefore also the number of bifurcation equations. 

In our considerations we distinguished the cases that the periodic orbit P has positive or negative Floquet 
multipliers, respectively. For positive Floquet multipliers we have discussed two different scenarios. First 
we have verified homoclinic snaking as it was previously observed numerically in our motivating example 
(1.2). Further we have described a nonsnaking scenario which has to our knowledge not yet been observed 
in systems in R 3 . In systems in R 4 however this effect was observed numerically. Though it is not clear 
whether or not in these examples this effect is due to the behaviour of the EtoP cycle as assumed in the 
present paper. For negative Floquet multipliers we confined to study the corresponding snaking scenario. 

For the detection of multi-around homoclinic orbit to E, these are orbits that follow the primary EtoP cycle 
several times before returning to E, also couplings near the equilibrium E have to be considered. At least 
if these investigations are combined with higher dimensional systems a Lin's method approach seems to be 
appropriate. However, numerical results in [25] give rise to the hope that those homoclinic orbit exist in our 
motivating example. 
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In Section 6 we verified numerically that system (1.2) satisfies the hypotheses which we used in our analysis. 
In particular we used a novel approach based on numerical continuation techniques to compute T. The 
computational results are very satisfying, the computed shape of T looks exactly like expected. This shows 
that Hypothesis 2.4(h) is true for our motivating system. 
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